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Abstract 

The newly developed single trajectory quadrature method is applied to a two- 
dimensional example. The results based on different versions of new perturbation 
expansion and the new Green's function deduced from this method are compared 
to each other, also compared to the result from the traditional perturbation theory. 
As the first application to higher-dimensional non-separable potential the obtained 
result further confirms the applicability and potential of this new method. 



PACS: ll.lO.Ef, 03.65. Ge 



1. Introduction 



Recently a new method has been developed in Refs.[l,2] to solve the low-lying quan- 
tum wave functions of Schroedinger equation using quadratures along a single trajectory. 
Based on the expansion on 1/g, where is a scale factor expressing the strength of the 
potential, Schroedinger equation can be cast into a series of first order partial differential 
equations, which is further reduced to a series of integrable first order ordinary differen- 
tial equations by single-trajectory quadratures. New perturbation series expansion and 
Green's functions of the wave equation are also derived based on this method, both for 
one-dimensional and N-dimensional cases. Some examples for one-dimensional problems 
have been illustrated in [1,2]. Recently this new method has successfully applied to solve 
quantum dot problem[3] and Yukawa potential[4]. 

In this paper, as the first complete illustration, this new method is applied to higher 
dimensional problems. Specially, it is shown in the paper how to derive the single trajec- 
tory based on Hamilton- Jacobi method. Schroedinger equation with a two dimensional 
non-separable potential is solved using this new method, based on two different versions of 
perturbation series, exponential and polynomial expansion, and using the Green's func- 
tion deduced from this method. The results show that these different versions are all 
equivalent. It is also shown in the paper that the result of the new method is the same as 
the one based on the traditional perturbation theory. It is much easier to obtain higher 
order perturbed wave functions using this new method, compared to the traditional per- 
turbation theory where the calculation of the second order perturbed wave function is 
already quite tedious. These results further confirm the potential of this new method in 
the future developments and applications. 

The paper is organized in the following way. In Section II. A brief review of the single 
trajectory quadrature method is given first. It is followed by a complete illustration to 
solve Schroedinger equation with a simple two-dimensional non-separable potential 



The same problem is solved using exponential and polynomial perturbation expansion, 
respectively, in Section III. In Section IV the Green's function method derived based on 
the new method is applied to the same example. Comparisons between different versions 
of this new method and to the traditional perturbation theory are given in Section V, 
together with a brief discussion and summary. 




(1.1) 
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2. Application to a two-dimensional example 

The method newly developed in refs.[l,2] provides a completely new way to solve the 
low-lying quantum wave function for one particle N-dimensional Schroedinger equation 
based on quadratures along a single trajectory. 

i7$(q)=E$(q), (2.1) 

where 

H^~W' + V{q) (2.2) 
is the Hamiltonian of a unit mass particle, and 

N Q2 

V' = E^- (2-3) 
The method consists three basic steps: 

1. For a potential V{q) > 0, a scale factor is introduced as 

V{c^) = gM<l)- (2-4) 

Expressing $(q) = e~^^^'^\ Both g'S(q) and the energy E are expanded in terms oll/g 
in the following way. 

^S(q) = <7So(q) + Si(q) + -S^ic^) + • • • (2.5) 

9 

E = gEQ + Ei + -E2 + ---. (2.6) 
9 

Substitute the expansion (2.5) and (2.6) into the Schroedinger equation (2.1) and equating 
the coefficients of (7^". After this step the second order partial differential Schroedinger 
equation changes to a series of first order partial differential equations of {Sj} and {Ei}. 

2. The lowest order equation of So in the series of equations, i. e., 

(VSo)' = 2v (2.7) 

is equivalent to a Hamilton-Jacobi equation in classical mechanics with —v as its potential 
and e = 0"^ as the total energy. It could be proved that the solution So can be written as 

T 

So(q) = / [^q^ - i-vmdt, (2.8) 

To 

where the integral is along the trajectory q(i) satisfying the classical equations of motion 

q(t) = Vv (2.9) 
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and the energy conservation 

\k{tf - = 0+ (2.10) 

The solution of (2.8) determines a single classical trajectory Sq. 

3. Introduce a new set of variables So and 

a = (ai(q),a2(q) ■ ■ ■ttAf-i(q)) (2.11) 

satisfying 

Va,-VSo = 0, j = l,2,---,iV-l. (2.12) 

Changing (gi . . . qN) (So, ai, • • • , ttAr-i), all {Sj(q)} become functions of Sq and a, i. e, 
{Sj(So, a)}. Based on (2.12) the series of first order partial differential equations of {Sj} 
reduces to a series of first order ordinary differential equations with variable So. These 
equations can be solved by quadratures along the single trajectory of constant a, i. e, 
along Sq. 

Now we turn to the example of two dimensional non-separable potential. To solve 
Schroedinger equation 

H^{x,y) = E^{x,y) (2.13) 

with 

H = -iV" + V(x, y) 

= \ , 

where V{x, y) = g'^v{x, y) and 

v{x, y) = ]^{x' + hY) + ^U{x, y). (2.15) 
we follow the three steps of the new method. 

1. Express $ = e^^^ and introduce the expansion (2.5) and (2.6) for gS and E. 
Substituting (2.5) and (2.6) into (2.13) the series of equations for {Sj} and {Ei} are 

obtained: 

(VSo)'= 2v, 
VSo-VSi= iV^So-^o, 

VSo-VS2= i[V2Si-(VSi)V^i, (2.16) 
VSo • VSs = i[V2S2 - 2(VSi) ■ (VS2)] - ^2, 

2. So can be expressed as 

T 

So(^, y) = j [2 + y') - (-^(^, vWt (2.17) 

To 
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and (x, y) satisfy the equations of motion: 



,2., _L "7du (2-18) 



In the following wc take U{x,y) = x'^y'^ as our example. To solve (2.18), further 
expansion of x and |/ on is introduced: 



X = Xo + jlXi + 11^X2 + 

y= I/O + /^yi + fj'^y2 + ■ 



(2.19) 



A series equations for different orders of /i is obtained: 



: Xo = Xo (2.20) 

yo = b'^yo 

jM" : Xi= Xi + 2xQyl (2-21) 
yi = b^yi + 2x^1/0 
: X2 = X2 + 2xiyg + 4xoyoyi 

y2 = &^y2 + 2xgyi + 4xoXiyo (2.22) 



We first solve eq. (2.20) and obtain 



yo = Cye^^ + dyC-^K 



(2.23) 



Introducing the initial condition: 

r x = 



at i ^ -oo (2.24) 



we have 



Xo = c^e^ 
yo = Cye^^. 



(2.25) 



The undetermined constants and Cy will be fixed after solving the series of equations 
of {xi} and {i/j}, by fixing the end point of the trajectory (x^, yx) at final time T. 

Now we solve equation (2.21) for (xi,yi). Substituting (2.25) into (2.21) we could 
obtain 

' - ^2^^+^)' , (2.26) 
Similarly, substituting (2.25) and (2.26) into (2.22) the solution could be obtained as 

- ^\ £i p(3+2b)t J 4 p(l+4b)tl 

•^2— 6+1 L 2(6+2) ^862(26+1)^ J (2 27) 

y^- £E£ir_£Lp(4+&)t J 4 J2+3b)t] ^' ' 

y2 — 6+iL8(fe+2)'^ ^ 26(6+1) (26+1)'^ J" 
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In principle, the whole series of {xi} , could be obtained by solving the series of 
equations in different orders of /x. In the following we will study the result up to the order 
of /x^, i.e., 

X = Xo + fxxi + fi^X2 ,^ 28) 

Now using the condition that (x, y) = {xt, Vt) si, t — T, we can fix and Cy up to /i^: 



„ _ r.r-^fl n I vl \ i j/?^(3b+l) i] 

L-x — •l^T'i I J- 26(1+6) 6(1+6)^^2+6^ 86(26+1) J J 

r — ■,,..^r~'>T ] ^ n ""t i ,,2 4 Ji/T i 4(3+^) ] 

Lj/ — (/TC <! J- 2(1+6) "^^ (1+6)^^26+1 ^ 8(2+6) J 



(2.29) 



It can be readily shown that the solution (2.28) satisfies the energy conservation condition, 
namely 

lii' + f) - ili^' + bV) + l^xY] = 0+. (2.30) 
Substituting (2.25)-(2.29) into (2.17) and taking Tq — — oo we obtain 

So = ^(X^ + + /L^^2 2 ^3,2 2^ ^2.31) 

"2^ 2(6+1) ^ 4(6+1)^ ^^6 + 2 26+1^ ^ ' 

3. Now we try to solve {Sj} and {Ei) using the single trajectory quadrature along Sq. 
After So is obtained, making the variable transformation, the series of equations (2.15) 
becomes a series of quadratures along the single trajectory So: 

Eq — |V^So|at 9=0) 
So 



Si(q) = Si(So,a)= / (^[iV% - £;o], 

El = |v2Si - (VSi)\, ,=o> 
S2(q) = S2(So,«) / -fL,{i[V2Si - (VSi)V^i}, 



(2.32) 



(VSo) 

£^2= |[V2S2-2(VSi)-(VS2)] 



at q=0' 



To perform the integration we first make the transformation 

So ,„ T 



dS 







J dt. (2.33) 



(^So)^ _ 

This can be derived by directly substituting (2.31) into the left hand side of (2.33). By 
direct integration over t, we obtain the following result of {Si} and {Ei} to the order of 
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//^ and 1/g : 

Eo= 1(1 + 6) 



■^1 ~ 4(l+b)'^'^ 6^ 4(1+6)^ U(2+6) "T" 6 (2+b)(l+26) 46(1+26) ^ 

G _ t^^ I y'^\ („2 I y'^\ (2.34) 



16(6+1)^ V-^ ^ 63^ 86(6+1)^ V-^ ^ 6 
8(1+6)^ I (l+26)(2+6) '^•^ ^ 6 ^ ^ 2 2+6 ^ 62(1+26) ^ / 

E2- -W^(6^ + 46 + l). 



Putting the above results together we have the approximate wave function and energy 
of the ground state, up to the order of /x^ and 1/g, as follows: 



gEo + Ei + g-^E2 



(2.35) 



The result shows clearly the two expansion series of /i and 1/g. In principle, this solution 
could reach higher accuracy by treating higher orders of /i and 1/g. This example gives us 
a general illustration of how to solve the ground state wave function for high-dimensional 
Schroedinger equation with non-separable potential. 



3. Two versions of the new perturbation expansion 
series 



For the same potential 



V{x, y) = g^[\{x' + 6V) + ^U{x, y)] (S.l) 



the problem can also be solved by a new perturbation expansion based on the new method. 
For this purpose, let us define 

g''liU{x,y) = eU{x,y), (3.2) 

i.e., e = g^^. Write 

V = VQ + eU (3.3) 

and 

Vo = \g\x' + hY)=9'v, (3.4) 

with 

T;o = ^(a;^ + 6V). (3.5) 
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The unperturbed Hamiltonian is 

^^0 = -^V^ + Vo. (3.6) 

It is very easy to solve 

HQe~^^° = EqC-^^'^. (3.7) 
The obtained trajectory is the same as (2.25), namely 

{ 

which satisfies the initial condition {x^y) = (0,0) at t = — oo. The constants and Cy 
are determined by the boundary condition {x, y) = {xt, yr) at t = T and given as 

( '^^'^'Z (3.80 



(3.7) gives 
For 

we introduce 



So = ^{x' + by'), Eo = ^{l + b). (3.9) 

= {Ho + eU)^ = (3.10) 

^ = e-^^'-^'x, (3.11) 

then expand E and x on Vfi*- For the energy i?, wc have the same expansion as (2.6). 
There are two ways to expand %, which will be shown in the following. 
1. Exponential expansion of x- 
Introducing 

X = e-»"'^^-^"'^3--. (3.12) 

Substituting (3.11)-(3.12) into (2.10), a series of equations of {Sj} and {Ei} are obtained 
as follows: 

(VSo)'= 2vo = x' + b'y' 
VSo-VSi= IV^So-Eo 

VSo-VS2= i[V%-(VSi)V^i + ef^ (3.13) 
VSo-VS3= |[V2S2-2(VSi)-(VS2)]-£;2 



This series of equations looks very similar to (2.15). However, there are two differences 
between them. First, the Sq here is the solution of (3.7), which is exact and very simple. 
Second, the equation for S2 in the order of has a factor eU due to the perturbation, 
while in (2.15) g^nU is lincluded in the potential v in the equation of Sq in the order 
of g"^. Now, we use the simple solution of Hq, i.e., (3.9), to define the single trajectory 
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So = + by"^) and perform the same procedure for (3.13), as we have done in Section 
II. To reach the same accuracy we have to take the expansion to the orders of and 
g^^. In the calculation x and y in the integrand are expressed as (3.8) and the integration 
/o^° rvs°i^ always transforms to j'^^dt. For the simple example with U = x^y^, the 
obtained results are 
S 

Si= 

^3 = 4(lT6) (^^ + V) -^3 = ^ 

*4 = ~ 4(1+6)^^ y ^2Tb + 1+26^ (3.14) 
^ 4(1+6)^ V 46(1+26) (2+6) (1+26) 6 4(2+6)^ 

8(1+6)^ Kl+26)(2+6) (^ + T'' + 2(2+6 + 62(1+26)-'] 



So= i(x2 + 6|/2) 1(1 + 5) 



^2 - 2(m)^^^y^ 



7 ^2 I _ (^2 , r ^ 

86(1+6)^ ^ b) 16(1+6)2 ^.-^ ^ b^) 



Remembering e = g^ix the above results are the same as (2.34). 
If we introduce 

g\U = g'^iiU (3.15) 

instead of eU , the same procedure could be taken. Only the second and third equations 
in (3.13) should be changed, due to the change of the order in g in the perturbation 
potential: 

VSo-VSi = iV^So-^o + AC/ 

vSo-vs2 = |[v2Si-(vSi)2]-£;i ^"^-'^^ 

The perturbation potential \U now enters the equation of Si in g^ order. To reach 
the same accuracy we need here the expansion only to the order of and g~^. For 
U{x,y) = x'^y'^, the results are exactly the same as (3.14) up to the order of and 
5f~^(i.e., and g~^). 
2. Polynomial expansion. 

Now, we introduce the expansion 

X= l + ^-lXi+^7"'x2 + --- (3.17) 

and perform the same calculation for the perturbations eU and gXU, respectively. Sub- 
stituting (2.6), (3.11) and (3.17) into (3.10), comparing the coefficients of the same order 
of g~'^, the series of equations of So, Si, {xi} and {Ei} could be obtained. 
For the perturbation eU we have 

(VSo)' = 2vo 
VSo-VSi = |V%-£;o 
VSo-Vxi + eU = El (3.18) 

VSo • Vx2 + VSi ■ Vxi - iV V + e^^Xi = ^1X1 + ^2 

VSq • Vx3 + VSi ■ Vx2 - ^ V2x2 + e^7x2 = E1X2 + E2X1 + ES 
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Based on the first equation of So, tlie same trajectory of So = ^{x"^ + by"^) is obtained. 
All {xi} could be solved by quadratures along the single trajectory Sq, which is again 
transformed to the integration of j]^^ dt. For U = x^y"^, the final results up to the order 
of and are 

So= \{x^ + hy^), E, = \{l + h) 

Si = 0, El = 

e 2 2 

•Vn — ^2 2 2 [4+6 2 I 46+1 21 p _ e 

-^3 - 8(1+6)"^ y b+6^ + 6(26+1) y i' -^3-46 /g^^gN 

•v.— „ f 4+6 4 I 36 ^2.,2 i 5„2-,2 i 46+1 4] 

16(1+6)^ 1 2(2+6)-^ (l+26)(2+6)'^ ^ ^ fe-^ y ^ 262(26+1)^ / 

^5 = 16(1+6)2 + F'' + 86(1+6)2 + — j 

"^8(1+6)2 I (l+26)(2+6) + T-) + 2I2+6 + 62(1+26))/ 

It is easy to prove that the polynomial expansion is equivalent to the exponential one. 

When we expand (3.12) as x = 1 - g-^S2 - 9-^83 ... + ^-^'^Sa - g'^Ss )^ + • • • 

and compare to (3.17), we have 

Xi = -S2 
X2 = —S3 + ^Sj 



I 1 3 

X3 — — S4 + 2(S2S3 + S3S2) — ^S2 



(3.20) 



N ^ 

Xiv = E (-1) feT • ( E Si, • Si2 . . . SjJ 

fe=l ■ k 

1=1 

Substituting (3.20) into (3.18) we return to (3.13) for S2, S3 ■ • ■. In fact, when we expand 
g-ff 82-9 S3... according to (3.14), the obtained expressions are equivalent to (3.19) up 
to the order of and g~^. 

Similarly, taking the perturbation gXU, the expansion gives the following series of 
equations: 

(VSo)' = 2t;o 
VSo-VSi = IV^So-^o + At/ 
VSo-Vxi + i[V2Si-(VSi)2] = E, 
VSo-Vx2 + VSi-Vxi-|V2xi + |[V2Si-(VSi)2]xi = E,xi + E2 

vSo-Vx3 + vSi-Vx2-iv2x2 + i[v2Si-(vSi)2]x2 = £;iX2 + £^2x1 + ^3 

(3.21) 

Performing the same procedure for our example with U = x'^y'^, expanding to the order 
of and g~^, the obtained final results are equivalent to (3.19) for g\ = e. 
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4. The new Green's function method 

Based on the new method the Green's functions are introduced in ref[2] and apphed to 
some one-dimensional examples. Now we are going to use the new Green's function to 
solve our two dimensional example. Assume that we have a two dimensional unperturbed 
Hamiltonian 

//o = -^V2 + yo(x,y) (4.1) 
with Vo{x,y) = g'^Vo{x,y). The corresponding Schroedinger equation 

Hoe-3^° = y£^oe"^^° (4.2) 

can be solved and gives the unperturbed ground state wave function e~^^° and energy 
gEo. Introducing a perturbed potential eU, the Schroedinger equation is now 

{Ho + eU)^ = {gEo + eA)^ (4.3) 

The perturbed wave function ^ = g-sSo^^ ^.^^^ |-,g obtained using the following Green's 
function expression: [2] 

^ = e-^^° + Ge{-U + A)^, (4.4) 



where 



G = e-^^oC(l + TC)-^e^^° (4.5) 



or 



X = 1 + C{1 + TC)-\{-U + A)x (4.5') 

and T = — |V^. To obtain the operator C in G we have to start from the trajectory 
So which is obtained by solving the unperturbed Schroedinger equation (4.2). After the 
variable transformation (x, y) — > (Sq, a) under the condition VSq • Va = 0, the operator 
C is defined as 

C = g-'9 hi (4.6) 

where 

(So|^|S;,) = ^(So-S;,) = | J (4.7) 

and 



^lo = [(VSo)T'- (4.8) 

When So is known the wave function ^ can be obtained according to (4.4) and (4.5) by 
quadratures along the single trajectory Sq. 
Now we turn to our example with 

Vo{x,y)^^{x' + bV)- (4.9) 
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As shown in Section III, the trajectory is given by (3.8)-(3.9) as 

1 



and 



So = -{x' + by') 



—T f 

X = xre e 
y = |/Te~*^e''*. 



Substituting the above trajectory into (4.6)-(4.8), we have 

1 



h 



2 



So ^2 + 

and ^ 

1 / dSo 



g J + by'^ 

Introducing the expressions of x{t) and y{t) in (4.11) into (4.13), we finally have 



1 } 
C = - / dt'. 

Q J 
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Considering the simple example of eU — ex'^y', we expand x as 

OO OO (X) CXD 

1=1 m=l 1=1 m=l 

Substituting (4.15) into (4.5'), we obtain the following equation: 

OO , OO OO OO 

E aix''+ E PrnV'"" + E E aimX'Y^ = 
1=1 m=l 1=1 m=l 

CeA-Ce(l + rC)~Vy2 



C OO 

+ Ce(l + TC)-^ E (Ax^' - x2('+i)y2) 
u=i 

OO , ^ 

+ E (A|/2--a;2y2(™+i))/5^ 

m=l 

+ E E (Ax^'y^m _ ^2il+l)y2im+l)^^\ 
1=1 m=l ' 

For any x^'y^™, after integration over t, we have 



Q J q2[1-\- mo] 



— OO 



and 

{-TC)x'Y'^ = ^2(7T^ ^^^^^^ ~ I)x2('-i)y2™ + 2m(2m - l)^''?/^^™-^)} . 
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Specially, we define, similarly to in ref.[2], 

C(-TC) V' = r|-x^('-") (n<l) 
Ci-TCYy-"-- = C;^;)y'^"^-") {n < m) ^^■'''> 

and 

C(-TC)'a;2' = CT]^^^^ 
C{-TCrjf"^ = CtL,(,) (4.20) 

The expressions of the F's are given in Appendix. 
Expanding 

(1 + TCy^ = 1 + (-TC) + {-TCf • • • , (4.21) 

applying (4.17) and (4.18) successively in (4.16), equating the coefficients of each term 
x^'y^"*, a series of equations for the unknown coefficients {«/}, {Pm} and {aim} are ob- 
tained. To solve these equations we need to further expand A and all the unknown 
coefficients on e: 

eA= eA(l) + e2A(2) + --- 

ai^ eai{l) + e'^ai{2) -\ 

Pm = e/?m(l) + e'/5™(2) + ■ • • 
aim= eairnil) + e^aim{2) -\ 

First the coefficient of the x^y^-term gives 



(4.22) 



A -rlj + E E (ArJ'^ - rl:,,^m+ihim 

1=1 m=l 

+ E(Ari(,)-r;;i,\)a, (4.23) 

oo , , 

+ E (Ar^,(,) - TllJ,^^)(3m = 

m=l 

Introducing the e-expansion, comparing the coefficient of the same order of e, we have, 
for e^-term, 

^(1) - = 4^ (4.24) 



A(2) = - E E (A(i)r,^i - r,^;i,J,+i)a,„(i) 

1=1 m=l 

-|(A(i)r.,,-r;;,:>,(i) 
- E (A(i)ri.,(,) - ry+i)/3».(i) 

m=l 

To obtain ^^(l), /^^(l) and a/^(l) we analyze the other terms in (4.16). After expansion 
on e, the -term gives 

oo , oo oo oo 

E«K1)^''+ E PmW^+E E a,^(l)xV"^ 

i m=l /=1 m=l (^4.^6j 

= -Cx2y2 _ C{-TC)xY- 
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Since 

_(^^2 2 _ 1 1 ^2 2 

^■^ y — 3 2(1+6)-^ y 
-C{-TC)xW = -^.^,{x^ + i) 



(4.27) 



we have 



ai(l)= -ij^iift, ai>i(l)=0 

-ii^mW' /5->i(l) = (4.28) 



4g2 6(1+6) ' 
■92(1+6) 

Substituting (4.28) into (4.23), we have 



aii(l)= -^2(iT6T' a/m(l)=0 for Z + m>2. 



Now we turn to the e^-term, keeping only the x^- and y^-terms we have 

ai(2)x2 + /?i(2)y2 = ^{1)C {ai{l)x^ + /?i(l)?/2) 

-C(-TC)'(ai(l)a;V + /9i(l)xV) 

+A(l)C(-TC)an(l)xV 

-C(-rC)3an(l)xV 



(4.30) 



Applying (4.17) and (4.18) successively, substituting (4.24), (4.28) into (4.30), comparing 
the coefficients of x^-term and y^-term, respectively, we have 



n, (')\ — 1 \ b+2 I 18 I 3 

"Iv^^ 1635 (1+6)2^ b ^ {2+b){\+2b\ ^ 

/^i(2) = TKnSTTZ^i^^ ~ 



16s5(l+6)2L 6 ' (2+6)(l+26) ' 2+6J ( A '^^\ 

1 r 2fe+l I 18 I 3 1 \^-^^) 

1635(1+6)^1- 63 6(2+6)(l+26) 62 (1+26) J" 



Following similar procedure for the other terms of the order e , we obtain the coefficients 
«ii(2) = i^4(T^[| + (i+26K2+6) ]> term; 
"2(2) = 3234(iS?2(2+6) ' x^-term; 
/32(2) - 32,.(i+6;W ' fory^-term; 
«2i(2) = g^3(i+y^(2+6) ^ for x^y2_terni; 

«i2(2) = 8,3,(i;t?(i+26) ' fo^ x^t/^-term. 



and 



"22(2) = o 2n^6^^ ' fo^ x^y^-term. 



?2(l+6)^' 

Comparing to (3.19) the above results are exactly the same up to the order of and g"^. 



5. Summary and conclusion 

The single trajectory quadrature method newly developed in Refs.[l,2] is applied to 
solve the ground state wave function and energy for Schroedinger equation with a two- 
dimensional non-separable potential. 
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Different versions of perturbation expansion, as well as the Green's function based on 
this new method arc tested using this example. The results are also compared to the one 
from the traditional perturbation theory (see Appendix B). The consistency of all the 
results proves the applicability and potential of the various versions of this new method. 
To solve a problem based on the new perturbation expansion method it seems better to 
use g\U as the perturbed potential. Compared to the potential ell, a faster convergence 
could be obtained; while compared to the potential g^^U, g\U could avoid the complex in 
calculating the trajectory Sq when /lU is included in v, as shown in section II. While the 
perturbation expansion seems to rely on the scale factor g, the Green's function method 
is more general and provides a much wider applicability. It is of interests to apply this 
new method to other areas of physics where Schroedinger equation with strong potential 
or with various kinds of perturbations needs to be solved. 

The authors would like to thank Professor T. D. Lee for his continuous instructions 
and advice. This work is partly supported by NNSFC (No. 19947001). 
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Appendix A 



In the calculation of (4.16) the following expressions are needed 

y^^2l _ 1 J_™2« 
g 21 

{-TC)x^' = 2^(2/ - l)a;2('-i) 
{-TC)y'"^= ^i(2m-l)y2(-i) 

rl - (2^-1)!! 

^ ':(^) ~ {2gy 
pi _ (2m-l)!! 
^ m,{y) - {2gbr 

Expression of some F^'^ applied in our paper are given in the following: 



pi,i _ _j_ 
^ 1,1 — 

r2^1 = 833(2+fe)(fe + 3) 
r^'^ — 1 {■ 3 _|_ 6\ 

^ 1,2 — 8s3(l+2b)U2 b) 

^2,2^ 32g*{l+b) [t+26(^ + ft) + 2+6(6 + 3)] 
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Appendix B 

For comparison the result for the same Hamiltonian based on the traditional perturbation 
theory is given in the following. 
For the Hamiltonian 



Ho = -Iv' + lg'ix' + bV) 



we introduce 



where 



H = Ho + eU 

eU = x^ + y^ 

V cu 



HmiO = (-1)™ e-«^^^^ 



and 



a; = g for s = a; 
5(6 for s = y. 

The unperturbed part gives 

where the unperturbed wave function and energy are 

E^mn = -{gm + gbn). 

Introducing matrix elements 

F^^n- jds ut\s) u^:\s), 
the values of F^^^, relevant to our problem, are given in the following table. 

m 

n 


2 

4 






2 


1 


1 


2ui 




1 


5 




2U! 





V3 
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Up to the order of e^, the perturbed energy and ground state wave function could be 
expressed as 

^i^olix, y) + e'^ip^olix, y). 

Based on the basic formula of the traditional perturbation theory, using the values of 
F^^^ given in the table, we have 

A£;« = <4°i|:rV|4°?,> 
- ^0 0^0 -4^2^,' 



(2) 




E ' 

m n 

-E 



^ (0) I 9 T| 



> 



p(o) _ p(0) 



■'0 



1 



(62 + 46 + 1), 



p(0) E^(O) 

(0), (62 + 6 + 1) 1 , 2 , 

- ^oo^^'^Hg^3i,2(i, + i) WTT)^'' 2^(6+1)^ 
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and 



4%M = E '{E 



/'p(O) p(0) wp(0) j^{0) . 
(0) |^2.,2L/,(0) ^,,(0) U2„2L/,(0) 



/p(0) p(0) \2 
l-C/Q ~ ^m' n') 

rpia) rpigb) jpig) rp(gb) rp{g) rpigb) T^igb) 

m m' n n' 00 00 m' On' 1 /_ .,\ 



n.' nl rn 



1 



16^2 63(6 + 1) 
1 

16^2 h\h^\) 
1 



32 63(6+1)2 



(262 + 85 + 1)^(0) (a;, ^) 

(62 + 86 + 2)4°)(^,?/) 
(562 + 346 + 5)V'f2(2^,y) 

(6 + 6)4°^^,^) 
(66 + l)V'^°)(x,y) 



(m's^ + n'gb)"^ 



16 62(6 + l)(6 + 2) 



16 63(6 + 1) (26 + 1) 

6 + 3 , (0) 



^^^•^^4 0(^,y) 

V3 36 + 1 (0) 



^fka;,!/)}. 



16 62(6+ 1)' 



Substituting the expression of iIj1^q{x, y), keeping the terms up to the order of g ^, writing 
the obtained wave function as 

oc e-^^(^'+^2/^) x{x, y) 
- e-^^°-x(a;, y). 



Comparing the coefficients of each x^y'^ term, the obtained x(a;, y) is exactly the same 
as the results based on the newly developed method shown in our paper. 
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